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Abstract 

Lack of any information on the CP violating phase 6c p weakens our abihty to determine neutrino 
mass hierarchy. Magic basehne of 7500 km was proposed to overcome this problem. However, to 
obtain large enough fluxes, at this very long baseline, one needs new techniques of generating 
high intensity neutrino beams. In this letter, we highlight the magical properties of a 2540 km 
baseline. At such a baseline, using a narrow band neutrino superbeam whose flux peaks around 
the energy 3.5 GeV, we can determine neutrino mass hierarchy independent of the CP phase. 
For sin^ 2^13 > 0.05, a very modest exposure of 10 Kiloton-years is sufficient to determine the 
hierarchy. For 0.02 < sin^ 2^13 < 0.05, an exposure of about 100 Kiloton-years is needed. 
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I. INTRODUCTION 



Neutrino experiments in the last decade have determined a number of neutrino parameters 
to a good accuracy. Among the currently unknown quantities are (i) the CHOOZ mixing 
angle 6'i3, (ii) the sign of atmospheric mass-squared difference and (iii) the CP violating 
phase 5cp- Experiments are being designed/constructed to measure these quantities. 

At present, there are three efforts to measure a non-zero value for 6*13 using reactor 
neutrinos as the source. In each of these experiments, the survival probability P{i'e ^e) 
will be measured using a pair of identical detectors, one close to the reactor and the other 
about a km away. The deficit in the far detector is a measure of sin^ 26'i3. Since these are all 
disapperance experiments, they should have very low systematic uncertainties, to measure 
the small value of sin^ 2^^i3. Double-CHOOZ will start taking data soon and it will see a 
positive signal if sin^ 2^13 > 0.04 1]. Daya Bay [2I and RENO 3] are expected to improve 
on this measurement. Daya Bay's final sensitivity extends up to sin^ 2^^i3 > 0.01 ^\. 

It is possible to determine sin^ 26*13 by measuring Pi^u^ — )• z/g) (P^e) in an accelerator 



experiment. T2K and NOi/A 6| experiments aim to do this. Even if these experiments 
see a positive signal, determination of sin^ 2^13 from their data will be subject to very 
large uncertainties because the probability P^e depends on all three unknowns mentioned 
above j^]- The reactor experiments, on the other hand, will give a clean measurement 
of sin^ 2^13 because P{i'e ^'e), at the relevant energies, depends only on this unknown 
neutrino parameter. 

We label the three neutrino mass eigenstates by their respective eigenvalues mi, 1712 
and m^. From the three masses, we can define two independent mass-squared differences 
A21 = m\ — m\ and A31 = m\ — m\. A21 is the mass-squared difference which drives the solar 
neutrino oscillations. It is known to be positive and its magnitude is much smaller than that 
of A31. Atmospheric neutrino oscillations are essentially driven by A31 whose magnitude 
is known but not the sign. If the neutrino masses follow the hierarchy 7713 ^ m2 > mi, 
called normal hierarchy (NH), A31 is positive. It is negative if the neutrino masses have 
the pattern m2 > mi ^ m3, called inverted hierarchy (IH). Present data allows both the 
possibilities. Determining the sign of A31 establishes the pattern (or hierarchy) of neutrino 
masses. In this letter, we propose a new scheme to realise this goal. 

Neutrino propagation through dense matter, leads to an effective mass-squared term 
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(usually called the matter term) in their Hamiltonian jsj. The interference of this term 
with the original mass-squared differences leads to the modification of the neutrino masses 
and mixing angles and hence their oscillation probabilities. The change induced by the 
interference, of course, depends on the sign of the original mass-squared difference and 
hence is different for NH and IH. By measuring this change in the oscillation probability 
P^e, induced by the matter term, we can determine the mass hierarchy. 

P^e depends on all three unknowns of neutrino physics. Hence it is very difficult to 
determine which of the three unknowns causes a change in P^e • Measurement of ^13 in reactor 
neutrino experiments leads to disentanglement of one parameter. The CP violating phase 
^cPi in principle, can be determined by measuring the difference in oscillation probabilities 
of neutrinos and anti- neutrinos. However, this method can not be used to disentangle 5c p 
and the matter term because the matter term changes sign under CP and induces a 'CP 
violation-like' change in the probabilities. Therefore a new strategy is needed to separate 
these two interlinked effects. 

A radical proposal was made sometime ago to disentangle 5c p from the matter term in 
P^e- The expression for P^e contains three terms and only one of them is dependent on 
5c p- If this term can be made to vanish, by an appropriate choice of neutrino energy and 
baseline, then it is possible to determine neutrino hierarchy without any information on 5cp- 
Calculation of the baseline, called magic baseline, gives an answer L ~ 7500 km [9] and it is 
independent of energy [lO|. Having sufficient neutrino fluxes at such a large distance from 
the source will be very difficult with the current accelerator technology. 

In this letter, we make an alternative proposal of a much shorter magical baseline. In the 
original magic baseline proposal, the condition that the 5c p dependent term vanish, holds 
both for NH and IH. This condition is quite restrictive and leads to such a large baseline. We 
propose an alternative condition which demands that the 5c p terms should vanish only for 
IH. This leads to a relation between the neutrino energy E and the baseline L. In addition, 
we also demand that P^e should be large for NH. This leads to a different condition on E 
and L. Solving these two equations gives us the solutions L = 2540 km and E = 3.3 GeV. 
At this energy, for this baseline, P^e is very small for IH and near the maximum for NH, for 
any value of 5c p- Thus a neutrino beam with its flux maximum at this energy can make a 
clean measurement of neutrino mass hiearchy, independent of 5c p- 
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II. CALCULATION 



A very good approximate expression for P^e, for three flavour oscillations including matter 
effects, is usually given as an expansion in the small parameter a = A21/A31. It can be 



written as 



sin2((l - A)A) 



+ a Ci 



sin((l-i)A) sin(iA) 



2 ^ sin^(AA) 



+ «^ C,^^, (1) 

where A = (I.27A31L/E) and A = A/A31. The matter term Q A (in eV^) = 0.76 x 
10~^p (gm/cc) E (GeV) eV^. p is the density of matter through which the neutrino prop- 
agates. Here A31 is given in units of eV^, L is in km and E is in GeV. The coefficients, Ci 
are given by 

Co = sin2^23sin2 2^i3 (2) 

Ci = cos 6^13 sin 26'i2 sin 26^13 sin 2^^23 cos(A + (^cp) (3) 

C2 = sin' 2^12 cos' ^23 (4) 

We note that only Ci among them depends on the phase Sep- 

P^e depends on all the three unknowns. Data on P^e from a single experiment leads to a 
degenerate set of solutions [12]. Data on P^e from experiments with different baselines can 



resolve some of these degeneracies 



13] . Here we assume that 6*13 will be measured in reactor 
neutrino experiments [iHsl, which resolves the degeneracies involving this parameter. Our 
proposal in this letter makes the hierarchy- ^c-p degeneracy irrelevant. 

In Eq. (Il]) the dependence on the matter term A is explicitly displayed. A31 is positive 
for NH and is negative for IH. A, on the other hand, is positive for neutrinos and is negative 
for anti-neutrinos. Thus, if we have only a neutrino beam, then A is positive for NH and 
is negative for IH. For anti-neutrino beam, the situation is reversed. In Eq. ([1]), the first 
term is the most sensitive to hierarchy but the second term provides a significant correction. 
As one varies 6cp in its range, the change in the second term can cancel the change in the 
first term caused by the change in hierarchy. In other words, there exist two degenerate 
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sets of solutions, (NH and 5cp = (V) and (IH and 5cp = ^2), both of which give the same 



value of P„e in a given experiment 



12| . This leads to hierarchy-^cp degeneracy and restricts 



our ability to determine the hierarchy. To overcome this, it was proposed to choose a 
baseline and energy for which sin(AA) = 0, so that the second and third terms vanish. The 
above constraint, for the first non-trivial zero, gives the magic baseline condition L ^ 7500 
km, independent of the energy j9|. The energy can be chosen by the condition that the 
oscillation probability be a maximum at this baseline. Therefore, one can indeed determine 
mass hierarchy at the magic baseline in a clean manner, independent of 5cp- However, one 
is now faced with the problem of obtaining large enough flux at this large a distance. 

We can make P^e independent of 5cp by choosing either sin(AA) = or sin((l — A) A) = 
[1^. The magic baseline uses the first condition which holds true for both NH and 
IH for the same L and is independent of E. The dependence of the second condition 
on L and E is different for NH and IH. We exploit this difference by choosing L and E 
such that sin((l — A)A) = holds for IH. This makes P^c for IH not only independent 
of 6c p but also very small because only the term in Eq. ([T]) survives. We impose the 
simultaneous demand that, for the same L and E, P^e for NH should be close to maximum. 
This leads to a substantial difference in P^g for NH and IH and enables us to determine 
the neutrino mass hierarchy for all values of 6c p, even for relatively small values of 613. 
The condition on P^c for IH translates into 1.27(|A3i| + A)L/E = vr whereas that for NH 
becomes 1.27(|A3i| — A)L/E = 7r/2. Solving the above two equations, we get L = 2540 km 
and E = 3.3 GeV, for IA31I = 2.5 x 10"^ eV^ [14|. Note that L is determined purely by A 
whereas E is determined by both A and A31. There is an uncertainty in L which is equal 
to the uncertainty in A. 

The large difference between NH and IH is illustrated in Fig. 1, where we have plotted P^e 
as a function of E for sin^ 2^13 = 0.02. The figure shows P^e as a band with 6cp spanning the 
entire 0° — 360° range, for both NH and IH. As claimed above, in the neighbourhood of E = 
3.3 GeV, P^e is about 0.002 for all 6cp for IH whereas it is 6 to 18 times larger (depending 
on 6cp) in the case of NH. This large difference can be measured if there is a substantial 
neutrino flux at E = 3.3 GeV. Note that the value of sin^ 2^^i3 chosen here is very small and 
is barely above the detectable limit of the experiments under construction. But even for 
such small values of 6*13, the configuration suggested here can make a distinction between 
NH and IH, independent of 6cp and using neutrino beam only. It is worth remarking that 
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FIG. 1: P^o as function of for L 2540 km and sin^ 2013 = 0.02. It is plotted for both NH and IH, in 
each case as a band. Within each band, 6c p varies in the range 0° — 360°. 



2540 km is the distance between the Brookhaven Laboratory and the Homestake Goldmine 



15|. The design for a neutrino superbeam, with flux peaking near 3.3 GeV, already exists. 



For example: NuMI beam with medium energy option at locations 7mr off-axis, has its peak 
flux at 3.5 GeV 

In this letter, we consider the following configuration. We assume that the z/^ source is 
located at Brookhaven which produces a NuMI-like beam with 0.8 megawatt beam power. 
This corresponds to 7.3 x 10^° protons on target (POT), with energy 120 GeV, per year 
j^. We assume a 300 Kton water Cerenkov detector at Homestake 2540 km away. We also 
assume that the orientation of the beamline is such that the detector location is 7mr off-axis. 
In our calculations, we take the NuMI fluxes and scale them appropriately to obtain the 



event numbers at this distance and for this off axis location 



16| . Our signal is electron 



appearance in the far detector, due to — > z/g oscillations. Single Tr*^ events produced by 
neutral current interactions form a potentially huge background to this signal. The visible 
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energy of the neutral current events is usually much smaller than the true energy of the 
neutrino and hence they can be suppressed by a large factor The remainder of these 
neutral current events, together with the electron events produced by beam v^. form the 



actual background. This 
of the unoscillated events 



Dackground was estimated in a previous study to be about 1% 



isl - We include this background in our study. We calculate the 
event rates for the energy range 1 — 10 GeV, in bins of width 0.4 GeV and smear the obtained 
event distribution with a Gaussian probability in energy with oe = 0.15E. We define the 
statistical between the event distribution for NH and that for IH, by 

Xstat N'^^ ' 

i=bins i 

Nf^ is the number of events for the true hierarchy plus the number of background events. 
Nf^ is the number of events for the test hierarchy, which is the opposite of the true hierarchy, 
plus the number of background events. The true hierarchy can be either NH or IH and we 
consider both possibilities. In calculating the event distributions, the following values of 
neutrino parameters are used: A21 = 8 x 10"^ eV^, IA31I = 2.5 x 10"^ eV^, sin^ 6^12 = 0.31 



and sin^^23 = 0.5 [l7[. We do the calculation for various different input values of sin^ 26*13 
starting from 0.1 and going down to 0.02 in steps of 0.01. 5cp is varied from to 360° in 
steps of 45°. 

We also assume a 2% systematic uncertainty in the neutrino flux and a similar uncertainty 
in detector systematics. The systematic uncertainty in the cross section is taken to be 10% 
{g]. These are taken into account through the method of pulls as described in 18-20|. In this 
method, the fluxes, cross sections etc are taken to be the central values in the computation 
of Nf^ but are allowed to deviate from their central values in the computation of Nf^. 
We assume that the k^^ input deviates from its central value by ak C,k, where ak is the 
uncertainty in this input. Then the value of Nj'^ with the changed inputs is given by 

npull 

Nr = Nristd) + E (6) 

k=l 

where Nf^{std) is the expected number of events in bin i for the test hierarchy, calculated 
using the central values of the fluxes, cross sections etc and npull is the number of inputs 
which have systematic uncertainties. The ^^'s are called the pu// variables and they determine 
the number of cr's by which the k^^ input deviates from its central value. In eq. (|6]), c\ is 
the change in Nf^ when the k^^ input is changed by [i.e. by 1 standard deviation). The 



uncertainties in the inputs are not very large. Therefore, in eq. ([6]) we consider only those 
changes in Nf^ which are linear in S,k- Thus we have a modified defined by 

X I4fcj = 2^ ^7tr + ?fc \<) 

i ^^i k=l 

where the additional term C,l is the penalty imposed for moving k^^ input away from its 
central value by cxk ^k- The with pulls, which includes the effects of all theoretical and 
systematic uncertainties, is obtained by minimizing x^i^k), given in eq. ([7]), with respect to 
all the pulls ^k'- 

(8) 

In addition to taking the systematic uncertainties into account, we have marginalized over 
IA31I, sin^26'i3, sin^ 2^23 and 6cp but held A21 and 612 fixed. In doing the marginalization, 
we assume that experimental uncertainties in IA31I and sin^ 26'23 are what they are expected 
to be from T2K (about 2%) [5]. It turns out that the marginalization over IA31I has on 
no effect but marginalization over 623 has a very significant effect, if IH happens to be the 
true hierarchy. We elucidate this point after discussing our results. 

The dominant term in P^e is proportional to sin^ 2^13, which makes the marginalization 
over ^13 the most crucial one in hierarchy determination. In the neighbourhood of oscillation 
maximum, the matter effects increase this term for NH and decrease it for IH, relative to 
its vacuum value. In shorter baseline experiments such as NOz/A, one expects to measure 
this increase/decrease and determine the hierarchy. However, it is possible to choose a 
6'^3, within the allowed range of 6*13, such that P^e(NH,6'i3) ~ P^c(III, 6*^3) [12]. In such a 
situation, marginalization over ^13 leads to a very small x^- In our proposal, the condition 
sin((l — A)A) = 0, makes P^e(IH) very small, in the neighbourhood of E = 3.3 GeV, 
independent of both 6c p and ^^13. In this energy range, P^e(NH) is close to oscillation 
maximum. It is also proportional to sin^ 2^^i3 and is quite large even for the very small value 
of sin^ 2^13, which is illustrated in Fig. 1. Therefore the oscillation pattern is very distinctive 
for each hierarchy and they can easily be distinguished for sin^ 2^^i3 > 0.02. We also note 
that the marginalization over 6c p ensures that this distinction between NH and IH exists 
even for the most unfavourable value of 6cp- Hence hierarchy determination is possible for 
the whole range of 6c p- 

Because of its importance, marginalization over sin^ 2^13 was done in different stages. At 
present, we only have an upper limit, sin^ 26'i3 ^ 0.1. Double Chooz jl| can measure this 
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parameter if it is > 0.04 and Daya Bay [2] can measure it for values > 0.01. Therefore, for 
input values sin^ 26*13 '■ 0.05 — 0.1 (that is if Double Chooz finds a positive result), we do 
marginalization only over this range. For input values sin^ 2^13 : 0.02 — 0.05 (expecting a 
positive result from Daya Bay but not from Double Chooz) we do marginalization only over 
this restricted range. 

It was mentioned above that the parameters A21 and 612 are kept fixed. These parameters 
occur only in the second and third terms of P^e- Given the smallness of a, the third term is 
very small and varying A21 and 612 within their ranges, changes this term by a very small 
amount. The second term undergoes much larger changes when sin 2^13 and 6cp are varied 
over their ranges. The change due to A21 and 6*12 is much smaller and can be neglected. 



III. RESULTS 



In our calculations, we assumed that only the neutrino beam is used. Neutrino beams 
have an advantage over the anti-neutrino beams because of the larger cross section and 
hence larger statistics. In table 1, we list the exposure, in Kiloton- years (Kt-yr), needed to 
obtain a 3a distinction between NH and IH if sin^ 26*13 > 0.05. In this case, one requires only 
minimal exposure (less than 10 Kt-yr) to distinguish the two hierarchies. This is independent 
of whether the true hierarchy is NH or IH. For 0.05 > sin^ 26*13 ^ 0.02, the results are shown 
in table 2. Here the needed exposure, for 3a distinction, sharply goes up with 6*13 from 7 
Kt-yr to 50 Kt-yr, if NH is the true hierarchy. If IH is the true hierarchy, the 3a exposure 
becomes about 92 Kt-yr. 

Note that the exposures, for NH being the true hierarchy, increase with decreasing values 
of sin^ 26*13 whereas, if IH is the true hierarchy, the exposures are more or less independent 
of the true value of 6*13. This feature occurs due to the marginalization over 6*13. The 
number of events in the region of peak flux of 3.5 GeV are strongly dependent on 6*13 for 
NH, whereas they are independent of ^13 for IH, as can be seen from eq. ([1]). If NH is 
the true hierarchy, N^^^ is computed using the input value of 6*13 but in computing N-^ 
we vary ^13 in its marginalizing range. Since the difference between these two numbers, 
decreases with decreasing 6*13, we need larger exposure to obtain the same xLin- If is the 
true hierarchy, then N^^ is computed using the input value of 6*13 and N^^ is computed 
with varying values of ^13. Around 3.5 GeV. N^^ is small and independent of ^13, whereas. 



9 



sin^20i3 (true) 


Exposure(NH) 


Exposure(IH) 


0.1 


2.93 


6.39 


0.09 


3.34 


6.04 


0.08 


3.94 


5.69 


0.07 


4.77 


5.38 


0.06 


6.04 


4.95 


0.05 


8.19 


4.55 



TABLE I: Exposure in Kiloton-years required for 3cr hierarchy discrimination in the case where 
both Double Chooz and Daya Bay see a positive signal. Second (third) column shows the results 
if NH (IH) is the true hierarchy. 



sin^20i3 (true) 


Exposure(NH) 


Exposure(IH) 


0.05 


7.32 


91.68 


0.04 


10.69 


86.80 


0.03 


18.97 


81.16 


0.02 


49.83 


71.50 



TABLE IL Exposure in Kiloton-years required for 3a hierarchy discrimination in the case where 
Daya Bay shows a positive signal, but Double Chooz does not. 

Nl^^ takes its smallest value for the smallest of the ^13 values in the marginalizing range. 
Therefore Xmin the exposure for 3a distinction, have a much weaker dependence on 
the input value of ^13, if IH is the true hierarchy. The small decrease in the exposure with 
decreasing values of 6*13 occurs due to the contribution of events in the energy range beyond 
4 GeV. For this range, the difference between N-^ and N.^^ increases as the input value of 
6*13 is decreased. Thus xLin increases and exposure decreases. 

For 0.02 < sin^ 26*13 < 0.05, the exposure is rather large if IH is the true hierarchy. This 
occurs due to the marginalization over 623. Without marginalization over this parameter, 
the required exposure is about 50 Kt-yr which is similar to the largest exposure needed if 
NH is the true hierarchy. With marginalization, the IH spectrum is computed with 623 equal 
to the input value of 7r/4 whereas the NH spectrum is computed with different values of 6*23 
in the allowed range jsl. Given that the number of events in the case of IH is very small, the 
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difference between the IH spectrum and NH spectrum is smaller, when the NH spectrum is 
computed with a smaller value of 023- Thus Xmin ^^^^ ^e smaller and we need larger exposure 
to obtain a 3a discrimination. 

The baseline 2540 km and the energy of peak flux 3.3 GcV were obtained based on the 
two conditions that P^e for IH should be independent of 6c p (which makes it very small) and 
for NH it should be close to maximum. This leads to a large difference between the expected 
number of signal events for NH and IH. However, the fluxes fall off as with distance. 
Therefore it is imperative to check that the baseline 2540 km is the optimum distance to 
determine mass hierarchy for a neutrino beam with peak flux around 3.3 GeV. We generate 
the event distributions for NH and IH for different baselines but with the same source (NuMI- 
like beam with medium energy option at 7mr off axis). The fluxes for different baselines 
are scaled by The input value of sin^ 26*13 is taken to be 0.02 and normal hierarchy 

is assumed to be the true hierarchy. The other parameter values and systematics are kept 
the same as in the earlier calculations. We compute the Xmin hierarchy discrimination 
for the baselines 1000, 1500, 2000, 2540 and 3000 km, assuming an exposure of 150 Kt-yrs. 
The results are shown in table 3. We note that the largest Xmm) larger by a factor of 3, 
occurs for 2540 km baseline, showing that this indeed is the optimum distance, for hierarchy 
discrimination, for the given neutrino flux profile. 



Baseline (in km) 




1000 


0.38 


1500 


0.54 


2000 


9.92 


2540 


26.20 


3000 


8.72 



TABLE 111: Baseline length L vs Xmin hierarchy discrimination for sin^ 2^13 = 0.02 if a NuMl- 
like beam with medium energy option at 7 mr off axis locations (peaking at 3.5 GeV) is used. 

We have highlighted some very interesting properties of a 2540 km baseline experiment 
and showed, through a simple numerical calculation, that such an experiment is well capable 
of determining neutrino mass hierarchy. In our calculations, we have included a conservative 
estimate of background, which we take to be 1% of the unoscillated events. Despite such 
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background, the setup we discussed is capable of hierarchy discrimination for even quite 
small values of ^^13. By imposing various kinematic cuts, the background can be suppressed 
at the cost of loss of some signal. This loss of signal can be compensated by having an 
increased exposure. However, the ability of the setup to determine the mass hierarchy will 
not be compromised because any such kinematic cut will lead to a larger signal to background 
ratio. 



IV. CONCLUSION 



In this letter we demonstrated the superior ability of a neutrino superbeam experiment 
with a baseline 2540 km, whose flux peaks in the energy range 3-4 GeV, to determine 
neutrino mass hierarchy. For sin^ 29is ^ 0.05, a very modest exposure of < 10 Kiloton- 
years is sufficient to distinguish the two hiearchies at 3a level. For 0.02 < sin^ 2^13 < 0.05, 
one needs an exposure < 100 Kiloton-years. These exposures are obtained analyzing the 
expected data from this superbeam set up only. If the data from this set up is analyzed in 
conjunction with the data from a reactor 613 measurement experiment, then the required 
exposures are likely to be much less. The set up we assumed is not hard to realize because 
2540 km is the distance from Brookhaven to Homestake and the technology for an accelerator 
beam, with peak flux in 3-4 GeV range, exists [6]. Such a set up, we believe, will have an 
excellent capability to measure not only small values of ^13 but 6cp as well. These issues 
are currently being studied. 
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